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ABSTRACT 


A collocation formulation is used as the basis for a unified approach to 
the various quasi-static aero-thermo-elastic problems. These problems 
include rigid and flexible load distributions, divergence, estimation of rigid 
and flexible static and dynamic stability derivatives, and the correction of 
wind tunnel data measured on flexible models. The formulation utilizes 


structural, thermal, and aerodynamic influence coefficients. 


The Aerospace IBM 7090 Computer Program No. LDOO03A provides the 
solution to the above problems. The program capacity is fifty collocation 


control points and ten values of dynamic pressure. 
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SYMBOLS 


Element of flexibility matrix 
Element of thermal influence coefficient matrix 
Reference semichord 


Element of oscillatory aerodynamic influence coefficient 
matrix 


Element of steady aerodynamic influence coefficient matrix 
Aerodynamic rolling moment coefficient 

Aerodynamic pitching moment coefficient 

Aerodynamic normal force coefficient 


Element of aerodynamic normal force coefficient distribution 
matrix 


Element of experimentally determined aerodynamic normal 
force coefficient distribution matrix 


Local chord 

Mean aerodynamic chord 

Average chord 

Local spanwise lift coefficient 

Element of force matrix 

Generalized reference deflection 

Reference deflection for initial deflection mode 
Element of detlection matrix 

Element of unit matrix 

Flexibility matrix normalizing constant 


Refezeace reduced frequency 


- iv= 


Rolling moment 

Element of mass matrix; bending moment 
Pitching patent 

Reference load factor 

Element of load factor matrix 

Dynamic pressure 

Divergence dynamic pressure 

Surface reference area 

Surface span measured from root to tip 
Reference temperature; torque 

Element of temperature distribution matrix 
Shear 

Element of aerodynamic weighting matrix 
Cartesian coordinates 

Coordinates of aerodynamic center 
Coordinate of pitching moment reference axis 
Normal force 

Denotes incremental value 

Angle of pitch 

Figenvalue 

Control point location as fraction of local chord 
Atmospheric density 


Frequency 


aV «= 


Subscripts 


H, T, N 


Oo 
r 


Matrices 


[ ] 
1} 
L | 
ie 
[| 


Aerodynamic 

Flexible 

As subscript to CZ, Cm, and C), denotes partial differ- 
entiation with respect to H, T, or N, respectively, e.g., 
CZy = 8Cz/dH; as subscript to Cz, merely denotes 


reference to H, T, or N 


Inertial; also used as dummy index, e.g., Ay; denotes 
width of ith strip 


Model; moment 
Initial value 


Rigid 


Square 
Column 
Row 
Inverse 


Diagonal 


Vie 


SECTION I 


FORMULATION OF PROBLEM 


A. Introduction 


The most general numerical formulation that can be made of the various 
aeroelastic problems is a collocation formulation expressed in terms of 
matrices of structural, inertial, and aerodynamic influence coefficients. 

Such a formulation of the flutter problem has been discussed in Ref. 1, and 
of the static aeroelastic problems in Ref. 2. The present report reconsiders 
the static aeroelastic problems by reviewing the derivations and extending 
the computational procedures of Ref. 2 in treating the problems of rigid and 
flexible load distributions, divergence, estimation of rigid and flexible static 
and dynamic stability derivatives, and correction of wind tunnel data meas- 


ured on flexible models. 


The basis for the present generalized aero-thermo-elastic analysis is 
found in recent advances in both aerodynamic and structural theories. The 
aerodynamic developments have proffered aerodynamic influence coefficients 
based ona number of theories, covering the entire Mach number range of 
current interest. The structural developments have provided structural 
and thermal influence coefficients for external and temperature loadings, 
respectively, for surfaces of arbitrary aspect ratio. The current state of 
the art of aerodynamic and structural influence coefficient analyses is sur- 


veyed in Refs. 3 and 4, respectively. 
B. Sign Convention 


We choose the NASA body stability axis system as our coordinate 
system: x forward, y starboard, and z down. Positive loads and deflections 
are in the same directions, and positive moments are given by applying the 


right-hand rule to the coordinate directions. 


= 


Cc. De rivation of Eq uations 


The force distribution {FI acting on a flexible lifting surface or body 
arises from two sources: the aerodynamic forces resulting from the attitude 
(as specified by the deflection mode {ht), and the inertial forces resulting 
from a distribution of load factor jn}. The aerodynamic forces are found 
from the weighted (the weighting matrix must be derived from experimental 
data on the force changes resulting from changes in attitude; see Refs. 5 and 6) 
steady aerodynamic influence coefficients (AICs), and/or a set of experi- 


mentally determined control point force coefficients 


{F,} = (qS/e)LW] [Ch] jh} + qs ice) (1) 


The use of the experimental force coefficients is generally necessary to 
account for the aeroelastic behavior arising from camber and/or wing-body 


interference. The inertial forces are found from the mass matrix 
= 2 


The total force distribution therefore becomes 


i) 


iF} 


iF} t {FI (3a) 


It 


(qS/e)[W] [Ch] {hf + qS icle)) +({M]{n} (3b) 


fy 


The deflection mode {h} is composed of the initial deflection mode of the 
rigid surface fn} the thermal distortion mode (which may be considered 
to be initially built in) which is found from the incremental temperature dis- 
tribution and the matrix of thermal influence coefficients (TICs), and the 


deformation mode of the flexible surface {hr} 
{hp= {hp tlapdfeh+ fhe}. (4) 


The deformation mode of the flexible surface is found from the total force 


distribution and the matrix of structural influence coefficients (SICs) 
{het =K{a]{F}, (5) 


where K is a normalizing factor to the SICs (introduced for convenience in 
studying variations in stiffness levels). Combining Eqs. (3b), (4), and (5) 


permits the solution for the deformation mode 
{hef= KA] [a] {Ff (6) 


where {F,} is the force distribution on the rigid system (again considering 


the thermal distortion to be initially built in). 
{rit = @S/ENWI(Cng] ({he} + [ag] {e}) + aS{Ck f+ iMd{n} (7) 


“A distinction is made later in the derivation [see discussion preceding 
Eq. (31)] between two sources of the initial deflection mode: that arising 
from any built-in twist or camber and that arising from the attitude (some- 
times called the additional deflection mode). The program computes the 
initial deflection mode as the sum of these, {nt = Ho{h,/Ho} + H{h/H}, and 
both sets of data are required. The distinction is not made at this point 
because it is not essential to the derivation. 
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where 
[A] = (1) - (KqS/eia](W] [Cy] )7’ (8) 


The total force distribution is found by returning to Eq. (3b) with Eqs. (4) 
and (6) 


{F}=(BI{F ys (9) 
where 


[B] = [1 + (KqS/ENW)[Cy,) CAI Lal (10) 


For the purpose of structural analysis, it is often convenient to convert the 
control point forces into the structural loads of shear, moment and torque at 
particular stations. A set of transformation matrices may be defined to 


compute 


tre shear 


{vp} =[V/F]{F} (11) 
the moment 

iM}=(M/F]{F (12) 
and the torque 

{T}=(T/F]{F} (13) 


.4- 


Each element in the load coefficient matrices is the load at a particular 


station due to a unit control point force. 


In certain longitudinal problems, it is desirable to adjust the attitude 
so that the flexible system sustains a specified total aerodynamic force. If 


we denote the specified total by Z, then 
LU({F} - [Mifn})= z (14) 


and if we extend the definition of {ht in Eq. (4) to include the change in pitch 


attitude, we have 
{hy} = {hy of - ae{ xt (15) 


The substitution of Eqs. (7), (9), and (15) into Eq. (14) leads to the solution 


for the incremental pitch angle 


Ae = (1/D\(C - Z) (16) 
where 
es LU([BFr} - [MI {n}) (17) 
D = (qS/e)I[B][w] [Chg] {x} (18) 
and 


{Fro} = (qas/eW] [Chg] (jn, 5} + [ap] {tt} + gs{cle)} +{MJ]{n} (19) 


Jes: 


The final force distribution in this case follows from Eqs. (7), (9), (15), 
and (19) 


{F}= [B] {Fr} - 40(qS/e)[B] [Ww] [Ch] { x} (20) 
The final deflection mode follows from Eqs. (4), (5), (15), and (20) 
{nfs {hy} - A0jx} + K[a}{F}+[ap]{t} (21) 


A variation of the foregoing is found in maintaining a constant total force 
equal to that given by the rigid condition (without the thermal distortion). In 


this case 
en = (e 
z= aSt|(A/ELW][Chglfrrot + tc) }) (22) 


and Eqs. (16), (17), (18), (20), and (21) provide the solution for the final force 


distribution and deflection mode. 


The possibility of divergence is seen in Eq. (6) as the inverse of the 
matrix [A] approaches singularity. The critical values of (Kqs/t) that yield 
a singularity are related o the eigenvalues of the matrix, and the divergence 


condition is specified by the eigenvalues of the equation 
Mhef=[al{w](Cy {he} (23) 
where 


\ = t/KSq ,., (24) 


from which the divergence dynamic pressure is found 


daiy = C/ MKS (25) 
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Any positive value of Vaiv shows a divergence possibility although 
usually only the first or second values have any practical significance. The 


corresponding eigenvectors of Eq. (23) are the divergence modes. 


The foregoing development of the force distribution on the flexible 
system provides the basis for estimating aerodynamic stability derivatives. 


The aerodynamic coefficients are defined by the following: 


C, = Z/qs (26a) 
= (1/aS)L1 ({F} - [MI {n}) (26b) 
c= Mast (27a) 
= -(1 /qS@.x - x9) ({F} - [M) {n}) (27) 
C, = &/qSs (28a) 
= (1/aSelyl ({F} - [MI{n}) (28b) 


The center of pressure is found from 


(xe - x )/c = -M|Z (29a) 
= -C_/C, (29b) 

y/s = 2L/Z (30a) 

= Cy [C, (30b) 


oy a 


grt A a 


OO ae rir eet ae eedity I qe, Sp, a C2 ee 


~ RES Tee 4 
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The force distribution for use in the above equations is found from Eqs. (7) 
and (9). To derive the stability derivatives, it is convenient to define a 
reference generalized deflection H, a reference temperature change T, and 

a reference load factor N. Itis further convenient to consider an initial 
deflection mode {hy} = Ho{hy/Ho} where Ho is a reference initial deflection as 
distinct from that caused by the attitude change (additional deflection) specified 
by H. * These modifications permit Eq. (7) to be rewritten in a more general 


form 
{F,}= as((ls f+ (1/E)[W][Cn,] (H,{h,/H} 
+ H{h/H} + Taq] jni))+ N[M] {n/N} (31) 


The initial coefficients and the stability derivatives are defined by 


C, =C, +C, H+, T+tC, N (32) 
Z Z.. Ley Zr are 
Ce IG. Gi MG 5 TC. ON (33) 
a EG H T N 
C,=C, +C, H+C, T+C, N (34) 
re) H T N 


where C7,,, Czy, and C7, are examples of aerodynamic, thermal and 
inertial derivatives, respectively. The substitutions of Eqs. (9) and (31) 
into Eqs. (26b), (27b), and (28b) and a comparison with the coefficients of 
H, T, and N in Eqs. (32) through (34) yields the following. The initial 


% 
See footnote, p. 3. 
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coefficients are 


Cc 


m 


“2. 7 LI}(B] iC, } 


= -(1/E)x - x JEBI {Cz} 


Oo 


where 


C= (1/8) Ly [B]{Cz,| 


(e 


{Cz f = {Cz f+ (Ho /ewilC, 1h ./H, } 


The aerodynamic stability derivatives are 


C 


m 


C, = (1/s)ly|[B] {Cc 
H 


where 


oP LTj[B] {Cz,,,} 


= -(1/c)Lx - x J[B] iC... 


zit 


{Corb = (/ENWI Ch h/H | 


The thermal derivatives are 


C2 =LII{BI{C,,} 


(35) 


(36) 


(37) 


(38) 


(39) 


(40) 


(41) 


(42) 


(43) 


EA EY TENT TH CRIS ACI ge NURI NRE tee irnsepaelie econ MMO! YN 


NCR IT EL, 


Con = -(1/E)Lx - x I[B] {Cz} (44) 
Cy = U/L yIBHC2} (45) 


where 
{Copf = (L/E0WILCy Jap] {t/T} (46) 


Finally, the inertial derivatives are 


a =L1) {Cz (47) 
Conn = -(1/Zjia - x 1 {Cznt (48) 
5 (/s)Ly {C244 (49) 
where 
{C, of = (K/eW) [Ch J(Al al MI {n/N} (50) 


The center of pressure for each of these loadings is found from Eqs. (29b) 


and (30b) using the appropriate coefficients or derivatives. 


The rigid aerodynamic stability derivatives may be estimated from 
Eqs. (39) through (42) by taking [B] = [IJ] (i-e., q =0). Equations (39) 
through (42) may also be used to estimate dynamic stability derivatives by 


using the (complex) oscillatory AICs which are defined by 


{F} = pwb@s[ W] [Ch] fh} (51) 
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By comparing, this detinition to the steady definition of Eq. (1), we see it is 
only necessary to replace [Cy,,] by 2k“ (€s/S)[C,] and to permit {h/H} to 
become complex in order to make all of the preceding development applicable 
to the oscillatory case. Since this is a significant feature of the present 
analysis, we digress to illustrate the calculation in some detail. We choose 
as an example the estimation of the symmetrical longitudinal stability deriv- 
atives. We begin with the general expressions for the lift and moment coef- 


ficients for transient longitudinal motion. 


C= Fe +C, (at/2v)+C, (9c /2V) 


Da ”q 
bees (az°/4v*) +, (oc? /4v-) + --- (52) 
D'a Dq 


C=C atC (at/2v)+c (@c/2V) 
a Da q 


: Cone (Se? /4V") " Si, OF /4v?) +--- (53) 


If we truncate the series at the points shown in Eqs. (52) and (53), and con- 


sider the case of harmonic pitching (6 = a) with amplitude a. then 


2 
C,=al1C, + ik(C +C -K“(c +C | (54) 
Z | Z ( ee a) ( Zn2, ea) 


2 
C -=a_IC + ikIC +C -k {[C +C ) (55) 
m | mo ( Mh, » ( mpe, MD, 


244% 


eb a at 


where the reduced frequency k = w/2V. In the case of harmonic plunging 
(98 = 0, a = h/V) with amplitude h 9’ then 


= na2/0[ ike - K-C a ik°C 2 | (56) 
a Da Da 


@ 2 =o 2 ‘a3 
C_ = no(2/0) ik, -k Cn - ik Cn 2 | (57) 
a Da Da 


The lift and moment coefficients can also be found from Eqs. (39), (40), and 
(42) by making the substitutions noted above (i.e., [B] =[1], and replacing 
[Ch] by 2k@(es/S)[C,]}). 


Cy = HU {C, (58) 
6 
Cat 7 HO/E)Lx ~ x9) {C24 (59) 
where 
2 
{C2,,t = 2k°(8/S)[W][C)] {h/H} (60) 


To appraise the various longitudinal derivatives it is necessary to evaluate 
Eqs. (58) through (60) for both motions, pitching and plunging. For pitching 
about station x). the reference deflection is taken as H = ao» and the addi- 
tional deflection mode is {h/H}= {x} - x}; for plunging, the reference deflec- 
tion is H = ho and the additional deflection mode is {h/H}= rf. The 
computational sequence is now evident. The evaluation of Eqs. (58) and (59) 
for the plunging motion yields two complex numbers which, when identified 
with the real and imaginary parts of Eqs. (56) and (57), lead to the four 
derivatives CZp,’ Cmp,’ CZp2,’ and Cmp2a (assuming, of course, that 
the static stability derivatives CZ, and Cm, have already been determined). 
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A similar evaluation of Eqs. (58) and (59) for the pitching motion yields 
another set of two complex numbers which, when identified with the real and 
imaginary parts of Eqs. (54) and (55), and with appropriate utilization of 

the angle of attack derivatives determined previously, lead to the remaining 
derivatives CZq Omg’ CZ p,q’ and C.,,_. Ina similar manner, many of 
the lateral-directional dynamic stability derivatives can be found from 


oscillatory AlCs. 


Our treatment of the stability derivatives provides a means of estimating 
a dimensionless spanwise load distribution and the spanwise variation cf the 
chordwise center of pressure on a lifting surface. This calculation requires 
the data on which the AICs are based, i.e., the location and number of the 
chordwise control points, and the widths of the strips into which the surface 
has been divided. The load distribution coeff. ients are needed to make these 


calculations. The flexible load distribution coefficients are given by 
CL f= (B]{C2f (61) 


where {c,} denotes any of the various distributions for the different loads 
except for the inertial loads; i.e., ic, } can be taken to represent any of 
{Cz} {Cory} or {Cz 7}; in the inertial case, {Con} is used directly, i.e., 
{Cz, f= {Cont The elements of {Cz + are listed in order for each surface 
strip as derived in the AICs and can be used to find the local dimensionless 
loading and center of pressure. Let us denote the sum of the load coefficients 


on the ith strip by 


ac, = Cc, (62) 
1 stripi 


and the moment about the leading edge by 


a. i » a7, : (63) 
1 stripi 
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where € is the dimensionless chordwise location of each control point aft of 


the leading edge. Then the spanwise loading on the ith strip is given by 


(e yc/%), = -ACy (s/ay;) (64) 


where Ay; is the load strip width. The local center of pressure is given by 


= AC [acy (65) 
1 1 


The problem of reducing wind tunnel data is ourlast consideration. The 
purpose of a wind tunnel test is the measurement of the so-called rigid loads 
on the configuration (unless the model is specifically designed for aeroelastic 
measurements). But since no model is completely rigid, some aeroelastic 
correction can be made to any wind tunnel measurement. We now consider 
this correction. The situation is described by Eq. (61), which may be 


written as 
{old} = [Bl {CL} (66) 


where icle)} is the set of force coefficients measured on the flexible model, 
[Bm] is the aeroelastic amplification matrix based on the model properties 
from Eq. (10), and {che} is the set of desired force coefficients on the rigid 


configuration. Obviously the rigid force coefficients follow from the inverse 
of Eq. (66) 


fol) = [Baal 7) {ols (67) 


The rigid force coefficients so determined then provide the starting point for 


the aeroelastic analysis of the prototype. 
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Fig. 1. Jet transport wing geometry. 
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SECTION II 


GENERAL DESCRIPTION OF INPUT 


A, Units 


All units are taken in the pound-inch system with two exceptions: the 
surface area is in square feet, and the dynamic pressure is in pounds per 


square foot. 


B. Example Problem 


We consider the static aeroelastic analysis of the five strip wing whose 
geometry is shownin Fig. 1. This is the jet transport wing analyzed by 
Bisplinghoff, Ashley, and Halfman throughout Ref. 7. We will use the follow- 
ing program options: trimmed loads with the structural loads suboption, di- 
vergence, aerodynamic stability derivatives (due to angle of attack), and 
inertial derivatives (due to symmetrical load factor). We proceed to assemble 
the data required as input to the program when it is to perform these options. 
The flexibility, mass, and aerodynamic influence coefficient matrices (all size 
10 X 10) are printed in the program output example (pp. 45 through 48) and 
need not be shown here. The flexibility and mass matrices are taken from 
Ref. 1; the aerodynamic influence coefficients are developed from the subsonic 
lifting surface method of Ref. 8 at a Mach number of zero. We shall list only 
the additional data needed for the representative options. A basic list of con- 


stants can be used for all of the options and is given below: 
= 225 inches 


1.0x 10-7 
564. 236 square feet 


N ® mM a 
Hl i 


il 


-41,919 pounds (this is the lift on one wing required in level flight 
neglecting any tail load or fuselage lift, and taking one-half the 
fuselage weight as 17,400 pounds) 


1,0 


an 
us 


As 


Pree NA, ee le 


OT TON ER FRI aOR len oR HRs = 


T =0 

N =1.0 

s = 500 inches 
xX, = 0 

A, = 0 


For the loads option we make the following assumptions: the aircraft is 
to be trimmed at three values of dynamic pressure: q = 400, 800, and 1200 
psf; the aerodynamic matrix needs no experimental correction (the weighting 
matrix is taken as unity); no lift has been derived experimentally (\cl)} = fof, 
size 10 X 1); there is no initial deflection mode (fn, /HG} = {ot, size 10 X 1]; 
the additional deflection mode corresponds to zero angle of attack ({n/es} = fo}, 
size 10 X 1); and the symmetrical load factor distribution matrix is a unit 
column ({n/N} = {r}, size 10 x 1). The control point chordwise coordinate 


matrix for use in the trim analysis is 


{x} = -63. 20 (inches) 


measured from the (unswept) 35 percent chord line as shown in Fig. 1. The 
final data necessary to compute the structural loads are the shear, moment, 
and torque coefficient matrices. Considering only the loads at the side of the 


fuselage, we have 
[V/F] = [I], size 1x 10, 
[M/F] = |45 45 141 141 223 223 323 323 413 413], 


[T/F] = [-20. 25 +81.00 -17. 85 +71, 40 -15. 80 +63. 20 -13. 30 +53. 20 -11..05 
+44. 20]. 
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The divergence option may be run simultaneously with the loads option 
since it requires only a pertion of the above information. We request two 


divergence modes and dynamic pressures to be found. 


The aerodynamic derivatives options will be carried out for the rigid 
case (q = 0) in addition to the previous list of dynamic pressures. The addi- 
tional deflection mode {h/H} used in the aerodynamic stability derivatives will 


be based on a one-radian angle of attack mode and is found from 


. 35343 
. 41371 
. 31154 
, 24616 
. 27576 Sper degree) 
. 10304 
. 23213 
. 92851 
. 19286 
. 77143 


{h/H} = -{x/57. 296} = 


The calculation of the load distribution and the rolling moment coefficient 
requires the following additional data shown in Fig. 1: the spanwise coordinate 


matrix 
Ly] =190 90 186 186 268 268 368 368 458 458] (inches), 


and the width of each strip and the dimensionless chcrdwise location of each 


control point on the strip. 


Strip No, Ay of ba 
l 93 0.25 0.75 
2 89 0.25 0.75 
3 91 0.25 0.75 
4 95 0.25 0.75 
5 87 0.25 0.75 
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It is possible to perform all or any selected number of program options 
with one input deck if all data are real, although, in the example problem 
taking {h/H} = {o} in the loads option does not permit calculation of the aero- 
dynamic stability derivatives and taking q = 0 in the stability derivative op- 
tion (to find derivatives for a rigid system) is not compatible with the trimmed 
loads option which requires q >0. To avoid computing irrelevant results, we 
use two data decks, reproducing the cards (data) in the first deck that can be 


used as a part of the second deck. 


C. Program Restrictions and Options 


1. Maximum matrix size is 50 x 50. 
2. Maximum number of dynamic pressures per deck is 10. 


3. Provision has been made to reserve a partition in the upper left- 
hand corner of the total AIC matrix for control points (and the flexibility and 
mass matrices must be compatible with this order of control points) whose 
aerodynamic forces may be neglected or found from an alternate theory to 
that used for the primary control points. This partition is termed the "external 
store'' region since external stores are an example of a source of additional 
control points requiring such special consideration. The maximum number of 


control points that can be reserved for external stores is 49. 


4, Maximura number of modes to be evaluated for the divergence op- 


tion is 25. 


5. For real and complex derivative calculations, the maximum number 


of strips is 20; the maximum number of control points per strip is 10. 


6. Multiple options may be performed with one deck except for the 


following restrictions: 


a. The structural loads option must be used as a suboption to the 


other loads options. 
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b. In any one input deck, all data must be for steady case (real) 


options or all data must be for oscillatory case (complex) options. 
c. The divergence option cannot be used in the complex case, 


7, Only one matrix of associated aerodynamic influence coefficients 
may be included in any one deck. When external stores are present, the 
matrix includes the AlICs for the stores; however, this area may be set to 


zero. 


8. In the oscillatory case, the following data must be input in complex 
form: AIC matrix, temperature mode, experimental forces, initial deflection 


mode, additional deflection modes, and load factor mode. 


9. Program options and significant equation numbers as assigned in 


Part C, Section I are listed below. 
a. Loads for constant root angle of attack; Eqs. (31) and (9). 
b. Loads for trimmed condition; Eqs. (31) and (20). 
c. Loads for constant lift coefficient; Eqs. (22), (31), and (20). 
d. Divergence, Eq. (23). 
e. Initial coefficients; Eqs. (35) through (37). 
f. Aerodynamic stability derivatives; Eqs. (39) through (41). 
Thermal derivatives; Eqs. (43) through (45). 
h. Inertial derivatives; Eqs. (47) through (49). 
i. Experimental data reduction; Eq. (67). 


j. Structural loads; Eqs. (11) through (13). 
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SECTION III 
_ DATA DECK SETUP 


A. Loading Order 
Data decks are loaded behind column binary deck LDOO3A. There are 


nineteen items that can be used inthe data deck. The deck is set up with the 


items entered in the following order: 


16. 


17. 


Heading card 

Data deck control card 
Constants cards 

Dynamic pressure card(s) 
Flexibility matrix 
Aerodynamic matrix 
Weighting matrix 

Thermal matrix 
Temperature mode 
Experimental force column 
Initial deflection mode 
Additional deflection mode 
Mass matrix 

Load factor column 
Chordwise coordinates 
Spanwise coordinates 


Strip widths and control point locations (percent chord) 
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18. Load coefficient matrices (shear, moment, torque) 
19. Second and successive ''additional deflection mode'' columns 


Items 1, 2, 3, and 6 must be present in all data decks; the presence (or 
absence) of the other items is determined by the option and data indicators 


used in the data deck control card (Item 2). 


B. Input Data Description 

1. The heading card may contain any alphanumeric characters desired 
in Columns 2 through 80. It is normally used for job title, engineer's name, 
and data. 

2. Data deck control card (FORMAT 1814): This card specifies which 
options are to be performed and also indicates the presence of certain data. 
The eighteen fields of this card are used as follows: 

Field 1. Column 4 contains a 1 if load option 1 (LOADS FOR 
CONSTANT ROOT ANGLE OF ATTACK) is to be 
performed; it contains a zero or blank if not to be 


performed. 


Field 2. Column 8 contains a 1 if load option 2 (LOADS FOR 
TRIMMED CONDITION) is to be performed, a zero 
or blank if not. 


Field 3. Column 12 contains a 1 if loadoption 3 (LOADS FOR 
CONSTANT LIFT COEFFICIENT) is to be performed, 


a zero or blank if not. 


Field 4. Column 16 contains the number of modes (25 
maximum) of the DIVERGENCE ogtion to be 


performed, a zero or blank if not. 


Field 5. Column 20 contains a 1 if the INITIAL COEFFICIENTS 


option is to be performed, a zero or blank if not. 


By A 


Field 6. 


Field 7. 


Field 8. 


Field 9. 


Field 10. 


Field 11. 


Field 12. 


Field 13. 


Field 14, 


Columns 23 and 24 contain the number of solutions 
for AERODYNAMIC STABILITY DERIVATIVES 
desired. This option permits varying the "additional 
deflection mode"! {n/H} when the other data are 
unchanged. The first mode is included as Data Item 
12, andthe second and successive modes are 


included as Data Item(s) 19. 


Column 28 contains a 1 if the problem is complex, a 


zero or blank if all data are real. 


Column 32 contains a 1 if the THERMAL STABILITY 
DERIVATIVES option is to be performed, a zero or 
blank if not. 


Column 36 contains a 1 if the INERTIAL DERIVA- 
TIVES option is to be performed, a zero or blank 


if not. 


Column 40 contains a 1 if a weighting matrix is 


included in the data deck, a zero or blank if excluded. 


Column 44 contains a 1 if a thermal matrix and 
associated temperature mode are included in the 


data deck, a zero or blank if excluded. 


Column 48 contains a 1 if an experimental force 
column matrix is included in the data deck, a zero or 
blank if exluded. 


Columns 51 and 52 contain the number of elements 
(control points) reserved for external stores, a zero 


or blank if none is reserved. 


Columns 55 and 56 contain the number of dynamic 
pressures (10 maximum) for which cycling options 


(repeated for each '"'q'') are to be performed. 
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Field 15. Column 60 contains a 1 if the EXPERIMENTAL 
DATA REDUCTION option is to be performed, a 


zero or blank if not. 


Field 16. Column 64 contains a 1! if the STRUCTURAL LOADS 


option is to be performed, a zero or blank if not. 


a 


Field 17. Columns 67 and 68 indicate the order of the system 


contained in the data deck (50 maximum). 


Field 18, Column 72 contains a | if the flexibility matrix is 


present in the data deck, a zero or blank if absent. 


In the subsequent description of input data we find it convenient to 
refer to above fields 1, 2, 3, 4, 5, 6, 8, 9, 15, and 16 as Options identified 
by the respective field numbers. The other fields contain control numbers 


and will be referred to as Item 2, Field (respective number). 


3. Constants cards (FORMAT 6E12.8): The two constants cards 


contain the items listed in the following order: 
a. CBAR: mean aerodynamic chord 


b. FLEXK: flexibility matrix normalizing constant (or scaling 


factor); when the matrix has not been normalized FLEXK = 1.0 
c. CAPS: surface reference area 
d. CAPZ: vertical force 
e. CAPH: generalized reference deflection 
f. CAPT: reference temperature change 
g. CAPN: reference load factor 
h. SMALS: surface span measured from root to tip 
i. CAPXO: coordinate pitching moment reference axis 


je CAPHO: reference deflection for initial deflection mode 
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4. Dynamic pressures (FORMAT 6E12.8): If six or fewer dynamic 
pressures are liste , onecard is used; seven through 10 pressures require 


two cards. The q card(s) are not required when performing Option 4 
(divergence) above, but unless Item 2, Field 14 is zero or blank, the q's 


must be input and must agree with the number in this field. 


5. The flexibility matrix [a] is input if its presence has been indicated 
in the data deck control card (Item 2, Field 18). If not included as input, a 
rigid case is assumed ([a] = [0]). The following card order and formats 


must be used when the matrix is input. 
a. Control card (FORMAT 1814) 


Field 1. Columns 3 and 4 contain m, the number of rows in 
the matrix (< 50). 

Field 2. This field is not used; it may be left blank. 

Field 3. Column 12 contains IFORM; if IFORM = I, the 
matrix elements are to be input using FORTRAN 
FORMAT 6E12, 8; if IFORM = 0, the elements will 
be input using column binary format. 

Field 4. Column 16 contains IROW; if IROW = 1, the matrix 
elements are input by row; if IROW = 0, the elements 


are input by column. 
b. Matrix elements (use format specified above) 


If IFORM = 1, and IROW = 1: use FORMAT 6E12Z. 8 and input 
the matrix element by row: each new row starts on a new card (line). 
If IFORM = 1, and IROW = 0: use FORMAT 6E12. 8 and input 
the elements by column with each new column beginning on a new card. 
If IFORM = 0, IROW must = 0: the. matrix elements are input 
using column binary format. This format should be used only if the data are 


available as punched-card output from appropriate computer programs. The 
elements must be punched by columns; Column | starts in Origin | and 
Column 2 starts in Location (1 +m). A TRA card must end the deck. (This 
transfer card has a 7 and a 9 punch in Column 1, Columns 2 through 72 are 
blank, and Columns 73 through 80 will contain the characters used for 
identifying and sequencing the deck. ) 
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6. The aerodynamic matrix [C, ] is always present in the data deck. 
It consists of two parts if external stores have been reserved in Item 2, 
Field 13, and it is real (steady case) unless indicated as complex (oscillatory 
case) inItem 2, Field 7. The [C,] matrix must be entered as follows with 
this exception: Items a and b (below) will be omitted when no external stores 


are reserved: 
a. Control card for store [C, ] (FORMAT 1814) 


Field 1. Columns 3 and 4 contain m, number of control 
points reserved for external stores (< 49); m = 0 
will direct the program to set the store [C, ] to 


zero. 
Field 2. This field is not used, it may be left blank. 
Field 3. Column 12 contain IFORM as defined in Item 5a. 
Field 4. Column 16 contains IROW as defined in Item 5a. 


b. Store matrix elements (format specified in Fields 3 and 4, 
above control card): Omit this input if m = 0 is used in the control card. 
When the external stores [c, ] matrix is input, use the format in the manner 
described in Item 5b. If the matrix elements are complex numbers 
(oscillatory case) consider the matrix size as m X 2m (imaginary parts of 
the elements form the even-numbered columns) and input as though all 


elements were real numbers. 


Cc. 1/k.. card (FORMAT 6E12,.8): The reference reduced velocity 
card is always present. For the steady case, it may be blank or contain any 
number, but for the oscillatory case it must be the 1/k,, associated with the 


aerodynamic matrix. 
d. Control card for surface (c,] matrix (FORMAT 1814) 


Field 1. Columns 3 and 4 contain m, the number of control 
points on the surface. (m = size of system minus 


number of store control point. ) 


-27= 


Field 2. Columns 7 and 8 contain L, the number of strips 
(partitions) in the surface. L= 1 for a full 


(unpartitioned) matrix. 
Field 3. Column 12 contains IFORM as previously defined. 
Field 4. Column 16 contains IROW as previously defined. 


e. The input listed below must be repeated for each partition 
(i= 1, L). 


(1) Partition size (FORMAT 18I4) 


Field 1. Columns 3 and 4 contain n, the number of control 
points on strip i (size of partition i). If L= 1, 


n= total number of control points on the surface. 


(2) Matrix elements in partition i (format specified in 


control card, Item 6d). 


When the matrix is for the steady case the elements in each partition are 
input in the same manner as the elements in the [a] matrix (Item 5b). In 
the oscillatory case, the elements are entered in the same manner as 


described for the stores oscillatory [c,] matrix (Item 6b). 


Note: Punched cards output from IBM programs based on the sources listed 
below may be used (with no alterations) for Items 6c, 6d, and 6e when the 
theory is appropriate to the problem. Punched cards output from the pro- 
gram based on slender-body theory may be used for Items 6a and 6b if the 
first card (1/k_) is removed from the deck. The list shown is complete only 


at this writing; additional compatible sources are currently under development. 


W. P. Rodden, E. F. Farkas, H. Malcom, and A. M. Kliszewski. 
''Aerodynamic Influence Coefficients from Incompressible Strip Theory: 
Analytical Development and Computational Procedure." Aerospace Corp. 
Report No. TDR-169(3230-11)TN-5, 3 September 1962. 
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W. P. Rodden, E. F. Farkas, H. Malcom, and A. M. Kliszewski. 
"Aerodynamic Influence Coefficients from Supersonic Strip Theory: 
Analytical Development and Computational Procedure.'' Aerospace Corp. 
Report No. TDR-169(3230-11)TN-1, 1 August 1962. 


W. P. Rodden, E. F. Farkas, H. Malcom, and A. M. Kliszewski. 
"Aerodynamic Influence Coefficients from Piston Theory: Analytical 
Development and Computational Procedure.'' Aerospace Corp. Report 
No. TDR-169(3230-11)TN-2, 15 August 1962. 


W. P. Rodden, E. F. Farkas, and G. Y. Takata. ''Aerodynamic 
Influence Coefficients from Slender-Body Theory: Analytical Development 
and Computational Procedure.'' Aerospace Corp. Report No. 

T DR-169(3230-11)TN-6, 31 October 1962. 


7. Weighting matrix: The weighting matrix [W] is input if its 
presence has been indicated in Item 2, Field 10. The matrix is entered 
with formats identical to those given for the [C.] matrix (Item 6a,b,c, d, e) 
with three exceptions: (1) the 1/k,. card (Item 6c) is omitted; (2) in 
repeating Item 6a if m = 0 the program will use a unit matrix for the store 
[W]; and (3) the matrix is always real. It is not required that [C, ] and [W] 
be input with the same format in the same data deck; i.e., one may be in 
column binary while the other uses FORTRAN format, and one may be 


partitioned while the other is a full matrix. 


8. Thermal influence coefficients matrix: The thermal matrix [a] 
is input when Item 2, Field 11 notes its presence. The matrix is required 
input for Option 8. The formats for entering [aq] are the same as those 


used for the flexibility matrix: Items 5a and 5b. 


9. Temperature distribution column (FORMAT 6E12.8): The temper- 
ature mode {t/T} always accompanies the above thermal matrix. For 
steady cases (real AICs), the {+/T} matrix is entered as one column, six 
elements per card. For oscillatory cases (complex AICs), the mode is 
complex with the imaginary part considered as a second column; tabulate by 


column with the second column starting a new card. 
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10. Experimental force column (FORMAT 6E12. 8): The experimental 
force column {cfs is input only if its presence has been indicated by Item 2, 
Field 12. The column (real or complex) is entered in the same manner as 


the column in Item 9. 


ll. Initial deflection mode (FORMAT 6E12.8): The initial deflection 
mode {h,/ HL} must be present to perform Option(s) 1, 2, 3, and/or 5. 
Again refer to Item 9 for entering the column (real or complex). If all 
elements in the column are zero, a sufficient number of blank cards may 


be used. 


12. Additional deflection mode (FORMAT 6E12. 8): The additional 
deflection {h/H} must be included for Option(s) 1, 2, 3, and/or 6. See 


Items 9 and 11 for inputing this column (real or complex). 


13. Mass matrix: The mass matrix [M] must be in the data deck 


when performing Option(s) 1, 2, 3, and/or 9. The matrix is input as follows: 


a. Control card 1 (FORMAT 18I4): Field 1 must contain NRC, 


order of the matrix (same as Item 2, Field 17). 


b. Control card(s) 2 through 7 (FORMAT 1814): Field 1 and 
consecutive fields in control card 2 and successive control cards contain in 


order LL., LH., LL., LH LL...., LH LL, is the row 


)’ l 2’ a ae NRC NRC* 
number where the first nonzero element appears in column i; LH. is the 
row number where the last nonzero element appears in column i. If only 

one nonzero element is in column i, its row number must be used for both 


LL. and LH.. 


c. Elements in [M] matrix (FORMAT 6E12.8): The elements are 


entered by column as shown below. Each column starts on a new line (card). 


Column 1: elements LL, through LH, 


Column 2: elements LL, through LH, 


Column NRC: elements LL... through LHypc 
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Any zero elements between LL, and LH. must be entered or their respective 
fields left blank. Ifallelements inacolumn are zero, then at least one 
element must be entered as zero (a blank card may be used); LL. and LH, 


for this zero element must appear in the above control cards. 


14. Load factor column (FORMAT 6E12.8): The load factor column 
{n/N} is included if [M] is used (Option(s) 1, 2, 3, and/or 9). If no 
correction is desired, input a unit column. The column is real for steady 


cases and complex for oscillatory cases. (See Item 9 for column input.) 


15. Chordwise coordinates (FORMAT 6E12.8): The control point 
coordinate matrix {x} or |x] is input for Option(s) 2, 3, 5, 6, 8, 9, and/or 
15. The matrix is entered only once (six elements per card) for any combi- 


nation of the above options. 


16. Spanwise coordinates (FORMAT 6E12.8): A spanwise control 
point coordinate matrix Ly} must be used with Option(s) 5, 6, 8, 9, and/or 


15. The elements in [ty] are the control point distances from the roll axis. 


17. Strip widths and control point locations (percent chord): The strip 
widths Ay. and the control point locations b.9J (given as a fraction of the 
local chords) are required for Option(s) 5, 6, 8, 9, and/or 15. These data 


have the following input format. 
a. Control card(s) (FORMAT 18I4) 


Field 1. This field must contain NSTRP = number of surface 
strips (maximum of 20 for above options). Do not 


include external stores as strips. 


Field 2. Field 2 and successive fields contain NCPT, = number 


of control points (< 10) on strip i. (i = 1, NSTRP) 


b. Ay. and Ej j= 1, NCPT.) (FORMAT 6E12.8): These data 
are entered by strips; each strip requires one line (card) when NCPT <5 and 
two cards when NCPT > 5. Start with the first card for strip 1 and input the 
data as follows: 


Field 1. Ay): width of strip 1. 
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Field 2. §,1: location of first control point for strip 1. 


Field 3. § 2: location of second control point for strip 1. 


Field (NCPT, + 1), § NCPT;: location of last control point 
for strip 1. (If NCPT> 5, then 6,6 would be in 
Field 1, second card). 


Repeat the above for each strip 


18. Transformation matrices for structural loads: The shear, moment, 


and torque matrices, [V/F], [M/F], and [T/F], are input only if Option 16 
is selected in conjunction with Option(s) 1, 2, or 3. The elements in any 
specific row of each matrix are the loads at a particular station due toa 


unit control point force. 


The matrices are input in the order in which they are listed (above 
pp.) with each using the same format (similar to the mass matrix but input 
by rows instead of columns). The format for [V/F] is given below; the same 


format is used for [M/F] and [T/F]. 


a. Control card 1 (FORMAT 18I4): Field 1 contains NVMT, the 
number of rows in [V/F]. Note the matrix size is (NVMT X NRC), where 
NRC agrees with Item 2, Field 17 and NVMT < 20. 

b. Control card(s) 2 through 7 (FORMAT 1814): Field 1 and 
consecutive fields of these cards contain LOW): LHIGH). LOW, LHIGH,, 
ee ae LOWyywr’ LHIGH wy wr’ LOW. is the column number of the first 
nonzero element in row i, and LHIGH, is the column number of the last 
nonzero element in row i. (i= 1, NVMT) 

c. Elements of [V/F] (FORMAT 6E12.8): The elements are 


entered by row inthe same manner as the columns in [M]. (Item 13c). 
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19. Second and successive "additional deflection mode'' columns 
(FORMAT 6E12. 8): Itis possibleto vary the additional deflection modal column 
{h/H} to obtain two or more solutions for real or complex aerodynamic 
stability derivatives. Use the input forrnat given for Item 9 (or 12) repeating 


the format for each new column. 


C. Table of Option Requirements 


To aid the program user in data deck setup, the tollowing table shows 
which of the above described 19 data items are required to make up a data 


deck for any option. 


DATA ITEM 

= 1 2 3 4 § 6 7 8 9 10 11 12 13 14 15 16 17 18 19 
1ijP P P P P PO OO 0 O P P P P A A A A 
2|P P P P P PO OO OO O P P P P P A A A 
2 31P P P P P PO 0 OO O P P P P P A A A 
je) 4i1P P P A P P O A A A A A A A A A A A A 
ls5|P PP POPOAAOPAAAPPPAA 
C.lp PP POPOAAAAPAAPPPAO 
8iP P P P O P O P P A A A A A P P P A A 
91/P P P P O P O A A A A A P P P P P A A 
15 j}P P P P P P O A A P A A A A P P P A A 
16 |<«——<— Set up for option (1, 2, or 3) being performed ——., A P A 

P = must be present; A = must be absent; O = optional inclusion; 


-= present when Option 16 is used with Option 1, 2, or 3. 


Note: The above table applies to options performed singly with the exception 
of Option 16 which must be used concurrently with Option 1, 2, or 3.. Where 
multiple options are to be performed with one data deck and an item must be 
present for one of the options and absent for another, that item must be in- 
cluded in the data deck. Also any data item indicated as present in the con- 
trol card (Item 2) must be included in the deck, even though it may not be used 


in the particular option(s) chosen. 
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D. Example Keypunch Forms 


Example keypunch forms are given on the following pages. Columns 73 
through 80 are reserved for data deck identification and sequencing. Only 
the cards with sequencing are used in the two data decks set up for the exam-~ 
ple problem; the lines with Columns 73 through 80 left blank are for descrip- 


tion of input. 
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SECTION IV 
PROGRAM OUTPUT 


Output for Each Option 


In addition to all input data, the following is printed for each option. 


Li 


Option 1 (for each dynamic pressure) 


a. 
b. 
C. 


d. 


Deformation mode 
Total deflection mode 
Final aerodynamic force distribution 


Total force distribution 


Option 2 (for each dynamic pressure) 


a. 


b. 


Cc. 


d. 


e. 


Incremental pitch angle 

Deformation mode 

Total deflection mode 

Final aerodynamic force distribution 


Total force distribution 


Option 3 (for each dynamic pressure) 


a. 


b. 


Aerodynamic lift 

Incremental pitch angle 

Deformation mode 

Total deflection mode 

Final aerodynamic force distribution 


Total force distribution 
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4. Option 4 (divergence) 


a. For each mode: mode number, eigenvalue, divergence 


pressure, and number of iterations 
b. Divergence modes (eigenvectors) 
c. Check eigenvalues 
d. Check eigenvectors 
5. Options 5, 6, 7, 8, and 9 (for each dynamic pressure) 
a. Distributed force coefficients 


b. Aerodynamic normal force, pitching moment, and rolling 


moment coefficients 
c. Total chordwise center of pressure 
d. Total spanwise center of pressure 
e. For each strip: spanwise loading and local center of pressure. 
6. Option 15 (for each dynamic pressure) 
a. Corrected experimental force coefficients 


b. Aerodynamic normal force, pitching moment, and 


rolling moment coefficients 
c. Total chordwise center of pressure 
d. Total spanwise center of pressure 
e. For each strip: spanwise loading and local center of pressure. 
7, Option 16 (for each dynamic pressure) 
a. All data listed under the associated option (1,2, or 3). 
b. Shear, moment, and torque at load stations. 
B. Example Printout 


The printout for the example problem is shown on th2 following pages. 
Some of the results shown may be compared with similar calculations presented 
in Ref. 7. 
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SECTION V 


PROCESSING INFORMATION 


A. Operation 


FORTRAN II MONITOR system 
B. Estimated Machine Time 


Because of the number of variables involved such as options, order of 
system, number of dynamic pressures, etc., it is not possible to give a 
practical method of estimating machine time. When the data are appropriate, 
considerable time may be saved by performing as many options as possible in 
One machine pass. Also, complete but independent data decks may be stacked 


to save program read-in time. 

C. Machine Components Used 
About 26,000 core storages 
Standard FORTRAN input tape (NTAPE2) 
Standard FORTRAN output print tape (NTAPE3) 


Four utility tapes (NTAPE4, NTAPE5, NTAPE6, NTAPE8) 
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SECTION VI 


PROGRAM NOTES 


A. Subroutines Used 


In addition to the main program (page 77) the following subprograms 


were used. 
l, MPRINT, prints matrix in matrix format (page 91). 


2. MMULTD, computes the product [c]} of the two matrices [a] and 
[B] (page 95). 


3. MREAD, reads a matrix in either column binary or FORTRAN 
format (page 98). 


4, BINRD, reads column binary cards (page 103). 

Bi. RDLN, reads and prints title card (page 113). 

6. MNVRSX, computes the inverse of a complex matrix (page 115). 
co INVERS, computes the inverse of a real matrix (page 118). 


8. SWEEPX, computes true mode and sweeps it from the related 


matrix (page 123). 


9. CENTER, computes the aerodynamic coefficients, centers of 


pressure, and spanwise loadings (page 128). 


10. LOADS, computes the three loads options and the structural loads 
suboption (page 133). 


11. DERIV, computes the flexible load coefficients used in the derivative 


options (page 143). 
12. MITER, matrix iteration (page 149). 
13. NPNRMX, vector normalization (page 160). 


14. DIVERG, divergent pressure calculations (page 165). 
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B. Generalized Tapes 


Input, print, and utility tapes in this coding are defined as logical units 
2,3,5,6,7, and 8; however, these may be altered by pi. .ug the desired units 
on symbolic cards HMO30827 through HMO30832 in the main program. 
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SECTION VII 
FLOW DIAGRAMS 
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‘ i 1 oe _é » 
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FLOW DIAGRAM — MAIN PROGR 


TAPE REWIND 
DIMENSIONS EQUIVALENCE ee FORMATS DEFINITIONS TAPES 


PRESENT 


2 


READ: YES 


CIOKR M,L, IFORM, IROW [ 


IFORM, IROW & 


NXST x NXST 


NO 


[w]}=[r] 


(NEXST x NEXST) 


store [A] 


NXSTx NXST 


OIne 


YES CALL 
DIVERG 


CALL 
THERMAL MREAD 
DATA M,N, IFORM, IROW A 
PRESENT ox [art] 


NO 
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GRAM - MAIN PROGRAM 


TITLI SREAL ; 
TAPE REWIND TITL22COMPLX CALL RDLN READ: READ: PRINT: 


DEFINITIONS TAPES MAXR=50 TI OPTIONS CARD CONSTANT CARD oo 
MAXQ=20 Mert CONSTANTS 


PRINT: 
READ: CALL MREAD store [a | FLEXIBILITY 
PRESENT M,N, IFORM, IROW [0] NTAPE 4 MATRIX 


[2] 


SiG: READ: READ. one. store [A | IF 
; NTAPE 4 WEIGHTING 


[a] 


NO 


YES 


PRINT 
WEIGHTING 
MATRIX 


[w] 


READ. READ: 
[w] IN BINARY [w] IN FORTRAN 


READ 
M,L, IFORM 


store [w][ch] 


NTAPE § 


CALL 
MMULTD 


{arrit} [or ]{t/T} 


YES 
READ: {c ie} 


-70- 
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FLOW DIAGRAM — MAIN PROGRAI 


INITIAL YES <ALL ADDITIONAL 
DEFLECTION MREAD CEFLECTION 
MOD {n / u} MODE 

hy/Ho 


SAS 


PRINT: 


P 
CHORDWISE ; READ: co! 

READ: COORDINATES be READ: NSTRP, NCPT; P 

Xj; izl,NRC Yj i=l, NRC ic1, NSTRP c 
NSTR 


REWIND 
NTAPE4, 


NTAPE5 


read [a] CALL MMULTO iF EXPER \vES 


NTAPES 


PRINT: 


CORRECTED CALL CENTER J td & ; REWIND (204) = 
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RAM — MAIN PROGRAM (CONTINUED) 


rr \ CALL MMULTD 
READ: 
MASS 
MATRIX MASS MATRIX {cnm} = 
NEEDED ELEMENTS 


{wass] {nu} 


REWIND READ: READ: 
TAPES NTAPE 5 [Ww] NTAPE 4 Ca) 
[¢H] 


PRINT: 
INVERSE 
ERROR 


a= [Aa] REAL 
B= [a] meas. 


PRINT: 
INVERSE 
ERROR 


CALL MNVRSX 


[ay 


REWIND iF IF DER. OPT 
LOADS OPT CALL LOADS CALL DERIV. 
oui Ny (sao wee 
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FLOW DIAGRAM —- LOADS SUBROUTINE 


— a> & ar 
COO oOz 


o 
é 


a. 


© [ 
{cia} {F:} * {} rs : mee = fu, sug} MeALTO 
i x {or} k ko” \enes di - 3 aes) Tag] \CH MENS 


o & 


e Ke 
o 2 ZZ 222+ , TURN OFF «© \\ | Tunnon F TURN ON 
)-fy3 a @aercis 
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SRAM — LOADS SUBROUTINE 


© © 


tc} = 
; woloes YET" h uted. wate 
not = sera {o} = {w] fc.) {c} _ 
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FLOW DIAGRAM —~- LOADS SUBROUT | 


© 


) 

00 Fiprafh aeee ON] SENSE. MULTO MUD inp 
SENSE MMUL L : 

ager +N {cnin} aN Ta] artes { po} a0) {F,} {rs} = [a] { o} A a 


mn (ox) 


MMULTO 


ell} f | + th}orty 


PRINT TITLE: 


PRINT TITLE: 
LOADS FOR 
TRIMMED 

CONDITION 


LIFT 
COEFFICIENT 


OF AT TACK 


i = (ost (ee ~ 
waco sane REWIND C sense °F 
p= [FMoN] {r} - ‘strenoul: {r} a <S> a & 
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2AGRAM — LOADS SUBROUTINE 
CALL 


ore 
hetated [AV fod 
G @ 
foetd | sae La sith 
“fort | LEO oot 


NT 


a73% 


CALL 
MMULTD 


{F}= (2) {r+} 


nd eof} 


© © 
(Ey, patie 
‘ o:(8) {wexx} i. [1) (0) {wenx 


Che 9 | de oe] "ate 
{Fo} : {F}-w{cnin i Lisi 


TURN ON ON 
SENSE SENSE 
> 
LIGHT LIGHT 
2 3 


FLOW DIAGRAM —~— DERIVATIVES SUBI 


27 

DIMENSIONS REWIND TAPES Pra a TURN ON Tul 

EQUIVALENCE TURN OFF COEFFICIENTS SENSE id 

COMMON SENSE LIGHTS REQUESTED LIGHT L 

FORMATS NX20 

TURN ON PRINT TITLE: PRINT TITLE: PI 
SENSE LIGHT 2 AERODYNAMIC THERMAL 

{o} : { h yn} STABILITY STABILITY ol 

DERIVATIVES DERIVATIVES 


OFF 


TURN ON 
SENSE 


READ 
NTAPE 5 


a=[w] [cr] 


LIGHT 
| 


5 
6 
CALL 
NOP(6) + SENSE MREAD i 
= NOP(6)—1 LIGHT 2 c 
ON o={vu} 
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SRAM — DERIVATIVES SUBROUTINE 


TURN ON PRINT 
SENSE LIGHT | TITLE: 


AERO 
STABILITY 
DERIVATIVES 


INITIAL 
REQUESTED {o} *Ho {ny/ Ho} COEFFICIENTS 


OFF 


PRINT TITLE: 
THERMAL 
STABILITY 

DERIVATIVES 


PRINT TITLE: 
INERTIAL 
DERIVATIVES 


TURN ON 
SENSE 
LIGHT 
I 


CALL YES 


— Seo 


PRINT: 
ADDITIONAL 


DEFLECTION 
MODE 
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SECTION VIII 


SYMBOLIC LISTING 


Sorne of the program FORTRAN symbols are defined below. 


FORTRAN Symbols 
A(I,J), B(I, J), C(I, J) 


ATTIT(I,J) 


C1lOKR 
CAPH 
CAPHO 
CAPN 
CAPS 
CAPT 
CAPXO 
CAPZ 
CBAR 
CNIN(I, J) 


CZRE(I,J) 
DELY(J) 

F (I,J) 
FA(I, J) 


FLEXK 
FMOM(J, K) 
FR(I,J) 
HADD(I, J) 
HR(I, J) 


Definitions 
Elements of matrix ''working" arrays; also elements 
of [a], [C,], and [W] 


Element of column matrix representing the product 
[ay] {t/T} 
1/k,. 
H 
H 
re) 
N 


S (reference area) 


Element of column matrix representing the product 
[M] {n/N} 


Element of experimental force column matrix 
Ay; (width of strip j) 
Element of total force distribution column 


Element of final aerodynamic force distribution 
column 


K 

Element in load coefficient matrix [M/F] 
Element of rigid force distribution column 
Element of additional deflection mode column 


Element of initial deflection mode column 
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FORTRAN Symbols 
IFORM 


IROW 

L 

LH(K) 
LHIGH(J) 
LL(k) 
LOW(J) 


NCPT(J) 
NEXST 
NOP(1) 
NQ 

NRC 
NSTRP 
NVMT 


Q(T) 


SHEAR(J,K) 
SMALS 
TCCP 
TORQU(J,K) 
TSCP 

X(1) 

Y (1) 

Z1(1, J) 


Definitions 


l if matrix is in FORTRAN format; 
0 if matrix is in column binary format 


IFORM 


IROW = 1 if matrix is input by rows; 


0 if matrix is input by columns 


Wow 


Number of strips (partitions) in the surface 
aerodynamic matrix 


Row number of last nonzero element in column k of 
mass matrix 


Column number of last nonzero element in row j 
of [V/F], [M/F], or [T/F] 


Row number of first nonzero element in column k 
of mass matrix 


Column number of first nonzero element in row j 
of [(V/F], [M/F],or [T/F] 


Number of control points for strip j 

Number of external stores reserved 

Option i, Field i of control card (data Item 2) ) 
Number of dynamic pressures read in * 
Order of current case 

Number of surface strips used for derivatives options 

Number of rows in [V/F],[M/F], or [T/F] 

(ith) dynamic pressure 

Element in load coefficient matrix [V/F] 

8 

Chordwise center of pressure 

Element in load coefficient matrix [T/F] 

Spanwise center of pressure 

Element in chordwise coordinate matrix 

Element in spanwise coordinate matrix 


Location of jth control point for strip i 


The symbolic listing of the program is shown on the following pages. 
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“INT dd GNV SOYVI SINVISNOD GNV Guvd SNOTIdD 


f (8°8139 HL) 7 CHI 


(8T‘t=l1 


*({1}dON) 


*T S23dV1N 3dVL INdDNI OV34 
KIdu W1v9 


(TSESdVINS 23d VIN) 


( “ANTINOVONS 3SUaANT WOsd NUNL3y sOUNHAHBE 
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(8°9T3SS *xLZ HT) / B°9T3S 
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CTSOCOWH 
65989 CWH 
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*X9E OHT / SdIULS HL SEIT SxHE 
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G3dVih ONIM3 
S+IdVIN ONIMIN 06 
Q@ = R3dVIN 
Lo = 93dViN 
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